DNA Research

The hybrid grouper 'Yushuban' (Cromileptes altivelis [female] x Epinephelus fuscoguttatus [male]) has shown promising performance in aquaculture trials, yet no complete reference genome has been available for this hybrid. We generated two gapless, telomere-to-telomere (T2T) haplotype-resolved genome assemblies of a single individual by integrating short-read, PacBio HiFi, Oxford Nanopore ultralong-read and Hi-C data. The maternal (Yushu_Calti_1.0) and paternal (Yushu_Efusc_1.0) haplotypes are 1,077.33 Mb and 1,083.36 Mb in length, with N50 values of 45.46 Mb and 45.55 Mb, respectively. Across all platforms, over 99.4% of reads mapped to each assembly, and Merqury-based quality values reached 57.7413 and 59.2068 (PacBio). BUSCO recovered 99.78% and 99.75% of complete core genes. We predicted 27,248 and 27,566 protein-coding genes, of which 98.52% and 99.23% were functionally annotated, and annotated telomeric motifs and repetitive elements (46.12% and 47.18% of each genome). These assemblies fill a gap in grouper genomic resources and provide a foundation for genetic dissection, parental genome comparison and molecular breeding of this hybrid.

Common tansy (Tanacetum vulgare), a significant agricultural invasive species, threatens forage productivity, livestock carrying capacity, and plant biodiversity across North American agroecosystems. Despite being a close relative to the ornamental Asteraceae Chrysanthemum, T. vulgare genomic studies are severely lacking, but remain critical to the development of novel weed management tools and the investigation of plant invasion evolution. In this study, we performed the assembly and annotation of the diploid genome (2n=2x=18) of an invasive North American T. vulgare genotype using long-read PacBio HiFi and Iso-Seq sequencing technology. The assembled haplotype genome comprised 9 pseudomolecules with a total size of ~4.8 Gb and a scaffold N50 of 546 Mb. In addition, both the plastid and mitochondrial genome sequences were assembled. BUSCO analysis revealed a high overall genome completeness with 69,343 predicted protein-coding genes (AED<0.6: 98%). A high-quality North American invasive reference T. vulgare genome, complementing a recent one assembled for a native European accession, is critical for the investigation of plant invasion biology, genome evolution, and the future development of novel genetic biocontrol weed management technologies.

The Olympia oyster (Ostrea lurida) was once a major food source and a key component of coastal marine ecosystems and culture on the US West Coast. Overfishing and habitat contamination in the early to mid-20th century severely depleted populations of Olympia oysters, and the species is currently considered a species of concern. Although multiple restoration efforts are underway, progress understanding the species biology and recovery potential has been limited by the historical lack of genomic resources. Here, we report a chromosome-level assembly of the Olympia oyster. The assembly was built using Oxford Nanopore reads and Hi-C to obtain a highly contiguous 1.03 Gb-assembly consisting of 10 chromosomes. RNA from gills, mantle, adductor muscle, and a mix of brooding eggs and larvae were also sequenced as 150 nt long paired-end reads which were used to guide the genome annotation. Over 52 K genes were predicted in the genome, of which 34K were protein-coding. A BUSCO analysis on the assembly showed 98.9% of the mollusca_odb12 gene set were complete (with 2.1% duplicated, and 0.8% missing). The short RNA reads were also used to identify tissue-level gene expression signatures in adductor muscles, egg and brooding larvae, and mantle. Development-related genes were highly expressed in eggs and brooding larvae whereas many xenobiotic metabolism genes were present in mantle tissues, highlighting mantles role in the interaction of the oyster with external environments. This new genomic resource will greatly advance our understanding of this keystone mollusk and provide a critical tool for guiding restoration efforts and the development of sustainable commercial production.

Mitochondrial genomes of angiosperms are unusual due to their large size, low gene content, and high structural variation. In the palm family (Arecaceae), mitochondrial architecture is highly dynamic, characterized by expansive sizes, fragmentation, and extensive recombination activity. This biological complexity presents a significant bioinformatic challenge, especially since only three complete palm mitochondrial genomes are publicly available to date. Here, we present the complete mitochondrial genome sequence for the Amazonian palm Euterpe precatoria, assembled de novo using PacBio HiFi long-read sequencing data. The assembly required extensive manual curation and graph tracing to resolve a complex recombination topology. The final resolved sequence spans 1,268,308 base pairs and is supported by a mean read coverage depth estimated at 746X. Annotation revealed a high degree of structural complexity, including 177 total gene features corresponding to 76 unique genes. This gene set includes 41 unique protein-coding genes (56 features), 32 unique tRNA genes (118 features), and three unique rRNA genes, reflecting the characteristic gene duplication and fragmentation documented in Arecaceae. The complete mitochondrial genome sequence is publicly available and provides a high-quality resource for understanding the mechanisms of plant mitochondrial genome expansion and structural evolution.

Macaranga tanarius (Euphorbiaceae), a myrmecophilic plant grown in subtropical countries, was identified as a source of Okinawan propolis, a honeybee product, which is glandular trichomes developed on the surface of the fruits. This tissue contains a high amount of characteristic geranylated flavonoids known as nymphaeols (nymphaeol A - C and isonymphaeol B), which exhibit strong antioxidant activity. In this study, we characterized the structure of glandular trichomes of M. tanarius, showing that they are sac-like structures containing drupelet-like aggregates of secretary cells, and a subcuticular cavity was also observed. From a cDNA library, we identified two prenyltransferase (PT) genes involved in biosynthesis of nymphaeols, and a PT gene grouped in a primary metabolic PT. Our biochemical experiments demonstrated that one of the former secondary metabolic PTs (MtPT1) is a B-ring-specific geranyltransferase for eriodictyol to yield nymphaeols B and isonymphaeol B. MtPT1 was unusual for a plant PT, an enzyme that yielded multiple reaction products. The other PT (MtPT3) transferred a geranyl moiety to the A-ring of eriodictyol to form nymphaeol A. Although dimethylallyl diphosphate (DMAPP) was not recognized as its prenyl donor substrate of eriodictyol, giving a B-ring geranylated eriodyctiol as the prenyl acceptor substrate, this enzyme used DMAPP to attach the prenyl moiety to the A-ring and produced nymphaeol C. These data suggest that M. tanarius has evolved its PT function to produce a diverse range of prenylated flavonoids using a limited set of genes.

Although Cinnamomum includes globally important spice and medicinal plants, genome-scale knowledge of Cinnamomum species native to Japan remains limited. The Nansei Islands are an approximately 1,200-km-long archipelago located between the Pacific Ocean and the East China Sea and represent a key biogeographic region shaped by deep-sea barriers such as the Tokara and Kerama Gaps. Here, we integrated chloroplast and nuclear genome data generated in this study with public datasets to examine the phylogenetic relationships of four Cinnamomum species native to the Nansei Islands and ten related species. Chloroplast genomes showed no major structural differences among the four island species, whereas nuclear genome data resolved two distinct evolutionary lineages in the Nansei Islands: C. yabunikkei, C. x durifruticeticola, and C. daphnoides formed a distinct Nansei Islands lineage separate from Sri Lankan, Indonesian, and Taiwanese lineages, whereas C. sieboldii was placed in a Taiwanese lineage. Ancestry estimation indicated that C. x durifruticeticola and C. daphnoides share a common ancestral component, while C. sieboldii exhibited a single homogeneous ancestral component, suggesting reduced ancestral diversity after divergence from a Taiwanese lineage. These findings reveal previously unrecognized genomic differentiation among Japanese Cinnamomum species and suggest that deep-sea barriers, dispersal, and long-term isolation shaped their diversification.

Thrips are known for being pests of many agricultural crops, as well as vectors of plant viruses. Yet despite their economic importance, genomic resources for this group remain limited, and no chromosome-level genome has been available for the genus Scirtothrips. One such species is the chilli thrips, Scirtothrips dorsalis Hood (Thysanoptera: Thripidae). This invasive and highly polyphagous pest is established across much of the world and transmits several tospoviruses. Genomic resources would provide a new means to understand the traits that make thrips difficult to manage. Here, we present a chromosome-level genome assembly for S. dorsalis. The assembly is 266.5 Mb, with 96.9% of the sequence anchored in 16 chromosome-scale scaffolds. It contains 20,557 genes, and repetitive elements account for 28.8% of the sequence, which is comparable to those of other sequenced Thripidae. A species tree based on 3,019 single copy orthologs placed S. dorsalis as sister to Scirtothrips hansoni, with Scirtothrips most closely related to Frankliniella, Megalurothrips, and Odontothrips. The mitochondrial genome is a single circular molecule of 15,207 bp in which the sequence of the published mini circle is integrated, in contrast to the bipartite mitogenome previously reported for the same cryptic species.

The mechanisation of agriculture has achieved major increases in labour productivity but has also led to an erosion of agrobiodiversity. Neglected and underutilized crops such as common buckwheat (Fagopyrum esculentum Moench) offer the potential to diversify food production systems but require major genetic improvement to compete with established crops. Although target traits for breeding such as more homogeneous ripening, reduced height, larger seeds and higher and more stable yields are well known, large-scale and field-based phenotyping data of diverse common buckwheat accessions across multiple environments are still scarce. Here, we tested the agronomic performance of a common buckwheat germplasm collection consisting of 170 accessions of diverse origin. The collection was tested at seven year-sites, and 19 agronomic traits related to phenology, growth habit and yield were assessed. Most traits showed high repeatability (0.31 < R2 < 0.98) and heritability (0.45 < H2 < 0.97), indicating genetic control of these traits. Coefficients of variation ranged between 6.19% and 46.15%, showing the high phenotypic diversity within the collection. By combining multiple traits, distinct ideotypes for different roles in the crop rotation were investigated, and accessions with promising agronomic profiles identified. Genome-wide association studies revealed marker-trait associations for 18 traits, allowing the identification of putative candidate genes influencing these traits such as a putative GRF-INTERACTING FACTOR associated with flowering time and biomass production. Our results show that the collection offers vast phenotypic diversity that can be exploited for the development of improved varieties tailored to meet specific agronomic purposes. The efficient utilization of these genetic resources should be supported by the application of genomics-based breeding approaches. The development of genetic markers and search for putative candidate genes mark important steps towards the establishment of efficient breeding schemes. Thus, these genomic resources provide a foundation for accelerating the genetic improvement of common buckwheat and developing agronomically competitive varieties.

Caraway (Carum carvi L., 2n = 2x = 20) is a medicinal and aromatic plant of the Apiaceae family valued for its essential oil with the major terpenoids carvone and limonene. While the essential oil content represents an important quality trait, low seed shattering and high thousand-grain weight (TGW) are key determinants of seed yield. As molecular tools remain scarce within Apiaceae breeding beyond carrot, we conducted a QTL mapping for these seed-related traits using two previously established F2 populations (Effective N(P1) = 236; N(P2) = 266). In P1 with the non-shattering maternal parent Rekord, QTL mapping revealed at least one shattering locus on LG03 (PVE: 9.9%). In addition, loci associated with essential oil (LG08, 10.5% PVE) and TGW (LG09, 11.1% PVE) were detected. In P2, we identified six essential-oil loci (LG02-06, LG10; combined PVE: 49.2%), with five loci affecting both carvone and limonene, and one (LG03) exclusively limonene. Additionally, a major QTL for TGW mapped to LG02 (PVE: 15.6%). For stalk attachment rate one minor QTL was detected (LG07, 8.5% PVE), while no QTL co-localized with the major shattering locus. To pave the path towards fine-mapping, we designed and validated 26 diagnostic PACE markers covering five key QTLs.

Bread wheat harbours low genetic diversity due to bottlenecks during speciation, domestication, and selection; wild relatives therefore provide critical reservoir for climate-resilient breeding. Aegilops mutica, one of the few self-incompatible diploids in the Aegilops-Triticum species complex, is a putative parent in the homoploid hybridisation event that formed the D-genome clade and is a proven donor of disease resistance, winter hardiness and grain zinc accumulation. Ae. mutica is globally endangered, with a reference-quality genome available for only a single Central Anatolian accession. Using HiFi sequencing, we generated genome assemblies from two accessions from opposite extremes of the species' range: western Turkiye and Armenia. The four haploid assemblies recovered ~93% complete BUSCOs and comprised ~81% repetitive sequence, predominantly LTR retrotransposons (~55%), closely matching the reference. Comparison against the reference identified 145,871-155,535 structural variants per haplotype, including an 8.4 Mb inversion on chromosome 6T. Between 96,860 and 105,472 structural variants segregated between haplotypes within individuals and classic domestication loci did not overlap with any large structural variants. These assemblies and structural variant catalogue extend genomic resources for Ae. mutica, providing a basis for targeted introgression into wheat and conservation of a species for which populations in the South Caucasus have severely declined.

Amaranthus tricolor is a climate-resilient orphan crop whose betalain pigments contribute to nutritional quality and stress tolerance. However, breeding efforts have traditionally relied on field yield trials due to the limited availability of genomic information and genetic resources. To address this gap, we first used the Illumina and PacBio sequencing platforms to assemble a novel reference genome of A. tricolor, using the Japanese cultivar Bayam. By combining this reference genome with a core collection developed in our previous study, we aimed to identify genomic regions associated with pigmentation-related traits across winter and summer seasonal environments. Betalains are antioxidant phytonutrients that contribute to abiotic stress tolerance, making the regulation of their biosynthesis under heat stress particularly relevant for developing climate-resilient crops such as A. tricolor. As a result, ATR1.0CH03_5131103, which is associated with ATR1.0ch03g000565, was identified under the winter seasonal environment. Meanwhile, transcriptomic profiles in this study revealed that ATR1.0ch03g000565 was differentially expressed under both control and heat stress conditions, suggesting that it may be associated with the regulation of CYP76AD1 expression under optimal conditions. Under the elevated-temperature treatment, transcript levels of CYP76AD1 and 5GT decreased in Am087, coinciding with lower betacyanin-dominated absorbance indices. cDOPA-like homologs were upregulated; however, pigmentation-related absorbance remained reduced. Further molecular experiments, such as CRISPR-Cas9 knockouts or overexpression studies, are necessary to validate the regulatory mechanisms of CYP76AD1 and 5GT under heat stress.

Phakopsora pachyrhizi, the causal agent of soybean rust, is among the most aggressive and economically impactful pathogens affecting soybean (Glycine max) production worldwide. As a dikaryotic rust fungus, P. pachyrhizi maintains its genetic information as two discrete haplotypes that are sequestered within separate nuclei for most of its life cycle as it reproduces asexually via the production of infectious urediniospores. Previous genomes of three P. pachyrhizi isolates collected in South America revealed a large genome with a high transposable element content, measuring up to 1.25 Gb and 93%, respectively. Although these efforts provided valuable foundational resources, the complexity of the genome resulted in highly fragmented assemblies. Here, we generated a highly contiguous, substantially improved, annotated genome resource from the P. pachyrhizi isolate LA04-1, collected during the initial outbreak of soybean rust in the continental United States in 2004. The assembly is haplotype-phased, with each haplotype assembled into 18 chromosomes totaling 572 Mb and 569 Mb in size. The isolate contains unique mating-type genes and does not share haplotype genomes with any available isolates with whole-genome sequencing datasets. This resource will be valuable to future investigations into the underlying mechanisms of virulence in this important plant pathogen.

The goat research community benefits from one of the most comprehensive livestock genome assemblies, ARS1. However, the rich diversity of goat breeds suggests that a single individual genome is not enough to represent the species full genomic content. Third-generation sequencing technologies have recently made it possible to generate very high quality assemblies that should enable to elucidate unique and/or complex phenotype determinisms depending on genetic structural variations. In the SeqOccin project, nine goat breeds were selected to represent French goat genetic diversity: Alpine, Boer, Creole, Fosses, Poitevine, Provencale, Pyreneenne, Rove and Saanen. High-molecular-weight DNA was extracted from blood samples and subjected to long read whole-genome sequencing. There reads were used to generate de novo assemblies and were aligned onto the public reference to call structural variants. PacBio HiFi and parental Illumina reads were used to assemble phased haplotypes of the Alpine breed. For the other breeds, PacBio raw (CLR) or ONT reads were generated. Our work demonstrates the feasibility of assembling genomes with a reasonable amount of data and appropriate software packages, delivering one high-quality assembly and 8 standard goat assemblies that pave the way for pangenome construction and structural variant functional studies.

In many flowering plants, self-incompatibility (SI) is a genetically controlled recognition system that prevents self-fertilization, thereby promoting outcrossing and maintaining genetic diversity within populations. The gametophytic SI system in grasses (Poaceae) is governed by two multiallelic loci, S and Z, each locus comprising of three genes governing the initial recognition between pollen and stigma. Identifying these genes establishes the molecular basis of SI in grasses, but how its unique two-locus architecture has shaped the evolution of the underlying loci and recognition genes remains largely unexplored. Here, we report targeted sequence capture of SI genes in 362 genotypes across nine forage grass species, revealing extensive allelic richness and sequence diversity. The SI genes exhibited elevated nucleotide diversity ({pi}) (0.134 < {pi} < 0.403) and great amino acid sequence divergence, ranging between 48.2-58.0% for the female determinants (sS and sZ) and 66.5-83.4% for the male determinants (SI-DUF247). The SI-DUF247 proteins also displayed trans-generic polymorphism, indicating that alleles predate speciation. Phylogenetic analysis of SI determinants showed congruence at both S- and Z-locus, consistent with co-evolution of recognition partners. Our results demonstrate that the unique two-locus SI system in grasses shares evolutionary characteristics with other SI systems, suggesting that their reproductive function is shaping their evolution. These findings also indicate that the grass SI loci are under long-term balancing selection to maintain recognition specificity while creating sufficient allelic diversity. Our study provides comprehensive sequence resources for SI genes and new insights into the evolution of the unique two-locus SI system of grasses.

Asian Passifloraceae occupy a distinctive evolutionary and ecological position within the family despite their relatively low species diversity compared with the Neotropics. Many Asian representatives are highly localized, small and often ephemeral climbers that persist seasonally through rhizomes, tuberous stems or dormant seeds, while considerable morphological plasticity may obscure species limits and evolutionary relationships. Xishuangbanna, southwestern China, lies at the convergence of Southeast Asian, Indian-Himalayan and southern Chinese floristic elements and represents an important regional centre of Passifloraceae diversity, with two species of Adenia and five species of Passiflora, including the recently described P. xishuangbannaensis and P. menghaiensis. We investigate the evolutionary relationships, plastome evolution and population genomics of these taxa using complete plastome assemblies, a concatenated dataset of 68 selected plastid protein-coding genes, 356 reconstructed nuclear Coding Sequences (CDS), and Genotyping-By-Sequencing (GBS) across the known wild distribution of P. xishuangbannaensis. Plastid, concatenated nuclear and nuclear coalescent phylogenies were strongly incongruent and differed from published classifications based on morphology and earlier phylogenies (ITS, ncpGS, ctyGS, ndhF and trnL-F). Gene-tree discordance and phylogenetic-network inference indicated widespread incomplete lineage sorting together with introgression. The two statistically best-supported network models placed P. menghaiensis, which showed high genetic similarity to the Indian P. napalensis, as a sister lineage involved in independent reticulation events associated with P. altebilobata, P. xishuangbannaensis and P. sumatrana. A pronounced biogeographic pattern was evident within Passiflora supersection Disemma, particularly in the P. siamica-P. cochinchinensis lineage across the Hua Line, whereas relationships among other Asian Passiflora blurred this boundary. This pattern suggests that diversification and reticulation within the Asian lineage may have overlapped temporally with the development of the biogeographic divide represented by the Hua Line. The concentration of early-diverging lineages around Xishuangbanna and the eastern Himalayan-Indian region further points to this area as an important centre in the early diversification of Asian Passiflora. Within this context, the highly localized P. xishuangbannaensis may represent either a relictual of an early Asian lineage or the surviving product of an ancient reticulate complex, contrasting with lineages such as P. siamica that subsequently expanded farther across Southeast Asia. Comparative plastome analyses revealed extensive structural and gene-content evolution within subgenus Decaloba, including pseudogenization, gene and intron losses, and repeated Inverted-Repeat (IR) expansions and contractions; substantial plastome modification also characterized supersection Disemma. Passiflora sumatrana exhibited a particularly distinctive IR expansion extending into the accD-pafII region, potentially characterizing this lineage. Adenia plastomes were generally more structurally conservative, although A. penangiana displayed an IR-boundary contraction involving the rps19-rpl2 region. Mitochondrial CDS and non-coding regions provided the best-resolved phylogenetic signal and revealed strong discordance with the chloroplast phylogeny, consistent with maternal mitochondrial inheritance and variable plastid transmission, potentially involving paternal or biparental inheritance and heteroplasmy. GBS analyses of P. xishuangbannaensis using complementary distant-reference, de novo and closely related-reference approaches showed that mapping to P. organensis provided the clearest SNP distribution and population resolution. ADMIXTURE, PCA and FST analyses revealed two strongly geographically structured populations with additional substructure and substantial differentiation (FST = 0.112). Linkage disequilibrium based demographic reconstruction using GONE2 indicated declining effective population size (Ne) toward the present in both populations, with a particularly abrupt recent decline in Pop1 and a more gradual decline in Pop2. Few candidate Self-Incompatibility (SI) genes could be confidently reconstructed, including limited recovery of highly variable docking regions. However, the intracellular region of an S-locus Receptor Kinases (SRK) contained a distinctive non-synonymous substitution predicted to alter protein configuration without major disruption of overall folding, potentially contributing to compatibility differences between populations. Together, these results reveal a complex evolutionary history combining historical dispersal, incomplete lineage sorting, reticulation, morphological plasticity, plastome structural evolution and strong population differentiation. The narrow distributions and recent demographic decline of these newly recognized Xishuangbanna taxa emphasize their conservation importance, while broader Asian taxon sampling and deeper genomic sequencing will be essential to resolve the evolutionary history of this discreet and poorly known Passifloraceae radiation.

Abstract: Picrotoxane-type sesquiterpene alkaloids, including dendrobine, are characteristic constituents of Dendrobium nobile, but the enzymes that produce their proposed skeletal precursors remain in completely characterized. Here, we combined tissue metabolite profiling with genome-guided screening and biochemical characterization of terpene synthases (TPSs). GC-MS analysis revealed variation in sesquiterpenoid profiles and dendrobine signals among tissues and developmental stages. Manual curation of the chromosome-level genome yielded 16 full-length TPS-a candidates. Functional screening in an FPP-supplying Escherichia coli system showed that three closely related genes on chromosome 1, DnTPS1~DnTPS3, produced copaborneol as their dominant product. The principal product of DnTPS1 was purified and identified by NMR spectroscopy, and a purified-enzyme assay confirmed its direct conversion of FPP to copaborneol. Fluorescence imaging in rice protoplasts and Nicotiana benthamiana leaves indicated predominantly cytosolic localization of DnTPS1. Guided by sequence comparisons and structural modeling, we tested 12 DnTPS1 variants and found that substitutions at W377 and I295 altered the relative distribution of sesquiterpene products. Phylogenetic and regional gene-order analyses placed DnTPS1~DnTPS3 in the same clade and identified related TPS candidates within corresponding genomic regions of six Dendrobium species. These findings establish that three D. nobile TPSs can form copaborneol and identify DnTPS1 residues that influence product selectivity, providing a basis for testing the proposed connection to dendrobine biosynthesis.

Parasitoid wasps develop in or on other arthropods and have evolved a wide range of strategies for exploiting their hosts. Ascogaster reticulata is a solitary egg to larval endoparasitoid wasp belonging to the subfamily Cheloninae of Braconidae. We constructed a genome assembly using long-read sequence data. The assembled genome was 280.5 Mbp, consisting of 209 contigs. The N50 value and GC content ratio were 13.9 Mbp and 30.43%, respectively. The percentage of transposable elements and repetitive regions was 52.43%. The predicted number of genes is 27,261. Merqury revealed that the QV value and completeness of the genome were 61.79 and 99.3917%, respectively. BUSCO analysis detected approximately 99% of the core genes. These results suggest that the genome data were of sufficient quality for use as a reference genome. Therefore, these data can be utilized as a genomic resource for comparative studies within Braconidae and for future investigations into the biology and evolution of Cheloninae parasitoid wasps.

Ancestry-related genomic variation may contribute to phenotypic diversity. However, its cellular consequences remain difficult to evaluate. To address this, we developed a human induced pluripotent stem cell (hiPSC)-based framework for investigating such phenotypes. We used the Japanese population as a model of Jomon-related and continental East Asian ancestry admixture. Jomon ancestry, estimated using the Jomon allele score (JAS), was associated with clinical measures, including renal-related traits and serum triglycerides in 34,152 Tohoku Medical Megabank participants. Peripheral blood mononuclear cell transcriptomics of donors from the upper and lower extremes of the JAS distribution (n = 10) suggested group-level differences in immune-related gene expression. To test whether the serum triglyceride association coincided with hepatic cellular phenotypes, we differentiated hiPSCs from these donors into hepatocyte-like cells. Neither intracellular triglycerides nor lipid-related gene expression differed between groups. These findings highlight the value of a genetically characterized hiPSC-based framework for investigating ancestry-associated phenotypes across cellular contexts.

The teleost whole-genome duplication (TGD) contributed to functional diversification of opsins. Some TGD paralogs, including those of parapinopsin (PP), Vertebrate Ancient (VA) opsin, and long wavelength-sensitive (LWS) opsin, show different absorption spectra and/or expression patterns. However, our knowledge of detailed evolutionary processes and mechanisms by which TGD contributed to opsin diversification is still limited. Here, we report that TGD paralogs of parietopsin (PT1 and PT2) are retained in extant species. While most species, including the zebrafish, have only PT1, several characins, including the red piranha (Pygocentrus nattereri), retain both PT1 and PT2. The Mexican tetra (Astyanax mexicanus) and catfishes have only PT2. To assess the degree of functional diversification between PT1 and PT2, we characterized spectral properties and expression patterns. Maximum absorption spectra differ slightly among PTs. Those of red piranha PT1, PT2, Mexican tetra PT2, and Japanese catfish (Silurus asotus) PT2 were located at 517 nm, 528 nm, 517 nm, and ~535 nm, respectively. Fluorescence in situ hybridization showed that (1) piranha PT1 and PT2 are expressed in the same pineal cells, and (2) they are also co-expressed with PP1. Based on these findings, functional diversification of opsins after gene duplication is discussed.

Premise:: Conifers comprise one-sixth of global forests and landscape genomic studies are urgently needed to understand their biology, yet relatively few conifer species have reference genomes. Targeted exon-capture is a practical alternative to whole-genome sequencing, but selecting loci without a reference genome remains a barrier. We developed a multi-tissue reference transcriptome for pi non pine (Pinus edulis) to ask how multi-tissue sampling improves probe design, and what is lost under a conventional single-tissue approach. Methods: We assembled a reference transcriptome from four diploid tissues (needle, phloem, root, vegetative bud) and one haploid tissue (megagametophyte), using the latter as a pseudo-guide. We designed a target-capture probe-set by prioritizing highly expressed P. edulis transcripts and transcripts matching a conifer exon database. Results: We recovered 76% of conserved plant genes with few duplicated transcripts. Filtering yielded 15,664 transcripts totaling 31.4 Mb, ~0.1% of the genome. Tissue-specific expression was substantial: 27.4% of transcripts occurred exclusively in non-needle tissues. Functional annotation confirmed complementary contributions, including seed-desiccation and germination genes from megagametophytes, and defense- and stress-associated genes from roots and phloem. Discussion: We built a transcriptome and efficient probe set using a roadmap suitable for other non-model plants with large genomes; keys to success were prioritizing complementary tissues and using expression-based filtering.